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We experimentally study the dynamical response of weakly-interacting atomic matter waves to 
a periodically pulsed, disordered optical lattice potential consisting of two overlapping standing- 
waves of incommensurate spatial periodicity. For periodic driving with a single standing wave 
only, we observe behavior consistent with the kicked-rotor model, namely ballistic spreading of 
momentum wavepackets at quantum resonances, and dynamical localization otherwise. However, 
adding the second standing wave can greatly modify these two effects. In particular, we find that the 
addition of disorder destroys dynamical localization when the driving is off-resonant and suppresses 
delocalization for a resonant drive. Our findings directly relate to quantum localization phenomena 
in low-dimensional systems, and illustrate the role of decoherence in the dynamics of driven quantum 
systems. 

PACS numbers: 05.45.-a ; 05.45.Mt ; 67.85.Hj ; 61.43.-j 
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In isolated ultracold atomic systems, the quantum na- 
ture of matter can be made manifest in sometimes strik- 
ing and counterintuitive ways. One example of distinctly 
quantum behavior arises in the dynamics of quantum sys- 
tems whose classical counterparts are chaotic, i.e. so- 
called quantum chaotic systems, where destructive inter- 
ference can suppress the onset of chaos [1 . The ^-kicked 
rotor model [1 has found particular relevance in exper- 
imental studies of driven quantum particles, as realized 
with cold neutral atoms periodically kicked by trains of 
optical standing wave pulses [2 . Whereas regimes of 
fully chaotic behavior with diffusive growth of the mo- 
mentum variable can be found in the classical system, 
the quantum analog generically displays dynamical lo- 
calization in momentum space due to quantum interfer- 
ence [H El S] , which has been shown [5 to be analogous 
to one-dimensional (ID) Anderson localization [6 in con- 
figuration space. Studies of the contrasting dynamics of 
classical and quantum systems, and in particular those 
addressing how classical behavior can arise in the pres- 
ence of noise [T-TT and nonlinearities [TT-TE' , are central 
to the understanding and control of quantum decoher- 
ence [16 . 

Foundational experimental studies of quantum chaos 
in driven cold atomic vapors ^ |17j have demonstrated 
dynamical localization. Moreover, recent progress has 
even allowed for the investigation of analogs to 3D An- 
derson localization in these systems [TB] [19] , as well as 
of the infiuence of noise and decoherence on quantum 
dynamics [2QH22J . Other recent studies of ^-kicked ro- 
tor dynamics with atomic matter-waves have allowed for 
clear observations of resonant dynamics [23] [24] when 
the pulse period corresponds to a discrete energy reso- 
nance in the system. Such resonantly pulsed matter- wave 
systems should also allow for studies of quantum deco- 
herence, owing to their connection to quantum random 
walks ^25j . 

In this manuscript, we explore the effects of controlled 



disorder on the dynamics of a quantum kicked rotor sys- 
tem, as realized with ultracold atomic matter waves in 
a periodically pulsed optical lattice. We observe that 
effects linked to interference in the quantum kicked ro- 
tor model, namely linear momentum growth at quantum 
resonances and dynamical localization for off-resonant 
driving, are largely destroyed due to the addition of dis- 
order. These observations suggest a breakdown in the 
quantum dynamics of driven particles under the infiu- 
ence of noise, corresponding to a quantum-to-classical 
transition in open systems. Moreover, our observation of 
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FIG. 1. Atomic matter waves in a periodically-pulsed stand- 
ing wave field, (a) A cloud of ultracold atoms is exposed to 
two incommensurate optical standing waves with wavelengths 
Ai,2 and spatial frequencies /ci(2) = 27r/Ai(2), whose intensi- 
ties are repeatedly pulsed on for a short time r at a periodic 
spacing T. (b) Time-of-flight momentum spectra, theoreti- 
cal (left) and experimental (right), of matter waves after N 
kicks from off-resonantly pulsed optical standing waves of ab- 
solute depths Vi = V2 = ^OEr, with Er the energy of atoms 
recoiling with momentum hki. 



disorder-induced transport for off-resonant driving has 
connections to Anderson localization in low-dimensional 
systems, as well as to transport in quasicrystalline mate- 
rials [26]. 

Our system consists of a nearly-pure Bose-Einstein 
condensate (BEG) of (1.4 ± 0.4) x 10^ ^'^Rb atoms in 
the |F, ttif) = |2, — 2) hyperfine ground state. The atoms 
are fully supported against gravity by a compensating 
magnetic field gradient, and are loosely trapped in a 
crossed-beam optical dipole trap of wavelength 1064 nm. 
The trapping potential has a mean harmonic frequency 
'^//o/27r ~ 20 Hz, and the BEG has a chemical poten- 
tial jj^/h ^ 0.3 kHz. To realize the quantum kicked ro- 
tor model, the atoms are driven by pulsed optical lat- 
tices oriented along the 2:-axis, as depicted in Fig. fl](a). 
We perform simultaneous pulsing by two optical lattices 
of wavelengths Ai = 1064 nm and A2 = 782 nm, and 
corresponding wavevectors /ci(2) = 27r/Ai(2). The lat- 
tice depths are characterized as Vi(2) = 8^2) ^R (with 
Er = h^k\/2M the recoil energy of the first lattice and 
M the atomic mass), and are calibrated by Kapitza- 
Dirac diffraction [27]. Subsequently, we use short pulses 
of length r = 2 fis (Raman-Nath regime), spaced at a 
pulse period T that is varied to investigate either reso- 
nant or off-resonant dynamics. After applying a train of 
N lattice pulses to the BEG, we immediately release the 
atoms and allow them to freely evolve in time-of-fiight 
(TOP) for 16 ms before performing absorptive imaging 
of momentum distributions, as displayed in Fig. [l](b). 

The system subject to a train of N lattice pulses 
is approximately [28 described by the ID Hamilto- 
nian H = -{fi^ /2M)d^ /dz'^ + S{t)ER, with S{t) = 
[si cos^ihz) + S2 cos^{k2z)] ^J^^ U{t/r - jT/r), where 
H is a unit-step function. By approximating the pulses 
as ^-functions, the Hamiltonian can be recast as that of 
a ^-kicked rotor [2^, H' = p^/2 + i^(0) J2f=i ^{^ - J). 
with K{(j)) = Ki cos(^) + K2 cos{r](j)) and r] defined by 
the ratio fe/Zci ~ 1.36, and rescaled quantities t' = t/T^ 
H' = {4.klT'^/M)H, (^ = 2kiz, p = {2kiT/M)p. Here, 
the conjugate variables ^ and p obey the commutation 
relation [0, p] = ii^, with an effective Planck constant 
/^ = SErT/H = ST/Tr (recoil time Tr = H/Er). The 
so-called stochasticity parameters i^i(2) = 45i(2)tT/T^ 
serve to delineate regimes of classically regular and 
chaotic motion. Resonant driving by the first (second) 
lattice is achieved when k/Att {rfn/An) is equal to a ratio- 
nal number [or in the parlance of Talbot interference ^24^ , 
when T matches a rational multiple of the Talbot time 

Tt = rW2(rwV)]- 

We begin our studies by investigating the effects of dis- 
order on dynamical localization in momentum space [2], 
using a pulse period of T = 36 /is that corresponds to off- 
resonant driving for each of the two lattices. For kicking 
with only a single lattice [i^i(2) = 2.3(0)], we find that 
the population remains trapped in momentum space with 
increasing kick number N ^ such that in Fig. El (a) we ob- 



serve apart from fiuctuations no sustained growth in the 
per-particle energy e (in units of Er). However, we see 
in Fig. [2] (b) that when we switch to a strong incommen- 
surate lattice potential {Ki = K2 = 2.3) the dynamical 
localization is destroyed and the population spreads out 
in momentum space, leading to an increase in kinetic 
energy. Such a deviation from the dynamics of the stan- 
dard kicked rotor, namely a transition from dynamical 
localization to delocalization, has been predicted for cold 
atoms periodically pulsed with incommensurate standing 
waves [To] , and it was proposed that for irrational values 
of T] this growth can be unbounded. 

To compare our experimental results to theory, we di- 
rectly simulate the matter-wave dynamics due to evolu- 
tion under H^ considering a basis of plane- wave momen- 
tum states coupled by the two light fields of the form 
|77i,n), having momentum Pm,n/^ki = 2m + 27177. We 
start with all population in the zeroth momentum order, 
i.e. iV^o) = |0,0), and then describe the total evolution 
from one pulse to the next by a unitary "kick" opera- 
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FIG. 2. Dynamical evolution of matter waves under off- 
resonant kicking in the presence of disorder, (a) The atoms' 
per-particle energy £ as a function of the number of kicks 
N. The atoms are kicked by only a single lattice potential, 
with stochasticity parameter Ki = 2.3, resulting in dynam- 
ical localization or a lack of growth in e. Black points are 
experimental data, the blue solid line and red dashed lines 
are simulated quantum and classical trajectories, as detailed 
in the text, (b) When an incommensurate, disordering lattice 
is added (K2 = 2.3), the dynamical localization is destroyed, 
and growth in the energy e is observed. (c,d,e) The rate of en- 
ergy growth with kick number, As/ AN, as determined from 
linear fits to experimental data and theory points. Filled black 
points are experimental, open blue disks are theory, and error 
bars represent the standard error of the linear fits. Shown are 
the cases of (c) only the first lattice with i^2 = 0, (d) only 
the second lattice with Ki = 0, and (e) both lattices present 
with the same strength of kicking Ki — K2. Red points in 
(c,d) are classical growth rates. 




FIG. 3. Energy growth under simultaneous off-resonant driv- 
ing from two optical fields, (a) Experimentally measured 
change in energy e from N — 1 to AQ kicks as a function 
of the stochasticity parameters Ki and K2 , for a pulse period 
T = 36 /is that is off-resonant for each lattice, (b) Change in 
mean per-particle energy from numerical simulations. 



tor tj [29]. After a train of N pulses, the final atomic 
wavefunction is given by IV^at) = t/^IV^o)- From the sim- 
ulated momentum distributions, we calculate both the 
mean per-particle energy e (in units of Er) and the rms 
momentum width cjp (in units of hki) [30]. As shown 
in Fig. [2] (a,b), we find in general good agreement be- 
tween the experimental measurements and the theoreti- 
cal quantum trajectories, which rely on no free param- 
eters. To demonstrate that the relative lack of growth 
in the single lattice case is due to dynamical localiza- 
tion and not to classical localization (KAM barriers sep- 
arating classically chaotic regions [31 ), we also show in 
Fig. [2] (a) the classical energy growth as averaged over 
10^ trajectories [32], which is nearly linear with N. 

Given the experimental restriction to relatively low val- 
ues of A/", and the appearance of non-trivial oscillatory 
behavior as in Fig. |2] (b), we abstain from attempting 
to discern fine details of the energy growth (e.g. sub- 
diffusive vs. diffusive). However, to gain some insight 
into the apparent disorder-driven crossover from localiza- 
tion to delocalization, we characterize the energy-growth 
per kick As/ AN using a simple linear fit to the energy 
evolution (for 1 < A/" < 60, allowing for an energy offset 
after the first kick) . The rates of growth as a function of 
the stochasticity parameters are plotted in Fig. [2] (c-e) 
for the cases Ki = 0, i^2 = 0, and Ki = K2. We see in 
Fig. [2] (c,d) no sustained growth for either of the single- 
lattice cases. We again rule out classical localization ef- 
fects by comparing to the linear growth rates of classical 
energy trajectories as plotted in Fig. [2] (c,d). When we 
kick simultaneously with both lattices, however, we ob- 
serve a transition to a regime of growth in Fig. [2] (e), 
both in experiment and in the simulations (with a region 
of weak transient growth near the transition, seen in the 



theoretical curves but not in experiment). 

To obtain a more complete picture of the dependence 
of dynamical localization on the two lattice strengths, we 
now more fully sample the available parameter-space of 
Ki and K2. In Fig. [3] (a,b) we show both experimental 
and theoretical results for the change in the kinetic en- 
ergy e between A" = 1 and 40 kicks, with good qualitative 
agreement. As expected, growth is minimal along either 
of the two axes, where only a single lattice is predom- 
inantly used, while significant growth is achieved when 
strongly driving the atoms with the two incommensurate 
lattices. A deeper understanding of the mechanism as- 
sociated with this growth may be found by considering 
analogies to atoms driven along one dimension by stand- 
ing waves at multiple temporal frequencies [H [181 EH ES] . 
In our scenario with multiple spatial frequencies, the 
kicking with only one lattice as in Fig.[2](c,d) would hence 
relate to Anderson localization in ID, while 2D localiza- 
tion is achieved by kicking equally with both lattices as 
in Fig. [2](e). Fig. |3](a,b) thus displays the crossover be- 
tween these two limits. Signatures of such a crossover in 
the Ki — K2 plane between dynamical localization and 
delocalization can also be obtained through an analysis 
of the kick operator tl that is independent of the kick 
number [32] , based on a stroboscopic Floquet state anal- 
ysis [5l [34]. Along the line Ki = 7^2, such an analy- 
sis reveals a transition to delocalized states at a value 
K12 ~ 2.2, in approximate agreement with our observa- 
tions in Fig. [21(e). Lastly, we note the formal similarity 
between our observation of delocalization in momentum- 
space induced by pulsed disorder and the recent obser- 
vation of disorder-enhanced transport in photonic qua- 
sicrystals [26], in which the increased density of states 
plays an important role in promoting transport. Simi- 
larly, when the atoms' momentum modes are coupled by 
both of the incommensurate optical lattices, the energy 
spacing between adjacent occupiable modes is greatly re- 
duced. 

We now turn to investigate the effect of disorder on 
the dynamics of a resonantly-kicked quantum rotor. For 
this purpose, the use of ultracold atoms with large de 
Broglie coherence length is vital, so as to avoid trivial 
decoherence due to spatial separation of the momentum 
components. To achieve a resonance condition for the 
first lattice, we set the pulse period to T = 124 /is, 
such that n/Aii ^ 1. Even for a weak optical potential 
{Ki = 1.6 ; si = 10), we observe in Fig. [4] (a) a linear 
increase in the momentum width of an ensemble of reso- 
nantly kicked atoms for up to A" = 20. We note roughly 
a factor of 2 difference in the growth rate from what we 
expect based on plane- wave simulations. This deviation 
can be explained through the finite size of the trapped 
atomic distribution (Thomas-Fermi radius ~ 9 /im. along 
z), which gives rise to an initial rms momentum spread 
Ap/hki ^ 0.03. We plot in Fig. [4] (a) the expected on- 
resonance momentum growth, as averaged over the initial 
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FIG. 4. Dynamical evolution of resonantly driven matter 
waves in the presence of disorder, (a) The momentum-width 
ap of the atomic distribution as a function of kick number N, 
for resonant driving of the first lattice at a pulse period of 
T — 124 /is. When kicked by only a single resonant stand- 
ing wave {Ki = 1.6), a near linear increase in Gp is observed. 
Black points are data from individual experimental realiza- 
tions, the solid blue line is a numerical calculation for an ini- 
tial plane- wave state, while the dashed red analytical curve 
takes into account the initial momentum- width of the atomic 
distribution due to finite-size. (b,c) Similar traces of the 
momentum- width as a function of N with added disorder, due 
to a strong incommensurate lattice potentials {K2 — 8,16). 
(d) The growth rate of the momentum- width, Aap/AN, as 
determined by a linear fit to the dependence on N, is observed 
to decrease with increasing strength of the incommensurate 
lattice K2. The growth rates determined from theory (open 
blue circles) are scaled by a factor of 1/2 to account for effects 
of the initial spread in momentum. 



momentum distribution of the sample ^ . To investigate 
the effects of disorder, we now add deep off-resonant, 
incommensurate lattices of strength K2 = 8 and 16 as 
shown in Fig. H] (b,c). While adding them leads to very 
large initial momentum widths after the first kick, the 
rates of momentum growth thereafter are found to be 
markedly less than in the single-lattice case. We again 
characterize the growth rate Aap/AN by a simple linear 
fit in the range 1 < A/" < 20, and these values are plot- 
ted in Fig. [4] (d) as a function of the second, disordering 
lattice strength K2. We observe a trend towards inhib- 
ited growth with an increasing amount of disorder, in 
agreement with the theoretical prediction. Such a break- 
down of coherent growth in the momentum width due to 
added disorder is consistent with the observed destruc- 
tion of dynamical localization, which also relied on quan- 
tum interferences. The observed sensitivity of resonant 
dynamics to spatial frequency noise mirrors the effects of 
decoherence by spontaneous emission on quantum reso- 
nances [2T] . 

In conclusion, we have studied quantum kicked rotor 
dynamics in the presence of controlled disorder, using an 
ensemble of ultracold atoms driven by two incommen- 
surate optical lattices. In the absence of disorder, we 
observe both a coherent growth in momentum when the 



pulse period corresponds to a quantum resonance, as well 
as dynamical localization in momentum space when the 
pulsing is off-resonant. We observe that the addition of 
disorder destroys dynamical localization and leads to mo- 
mentum growth for off-resonant driving, while also sup- 
pressing the growth in momentum space for resonantly- 
kicked rotors. Our observations highlight the fragile na- 
ture of effects relying on quantum interference, as well as 
the transition from quantum to classical behavior due to 
environmental coupling. 

Possible extensions of our study include the realization 
of kicked harmonic oscillators [36l El] with strong non- 
linear couplings, e. g. by studying the single-site orbital 
dynamics of a periodically-tilted Mott insulator with tun- 
able filling and interaction strength, or the study of the 
anomalous quantum Hall effect with a spinful version of 
the quantum kicked rotor [38], as may be realized with 
cold atom mixtures in a spin-dependent lattice [39] . 
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CLASSICAL DYNAMICS 

It is well known that classical localization can occur in kicked-rotor systems for small 
values of the stochasticity parameter, where the classical system is not globally chaotic but 
supports bounded orbits due to Kolmogorov-Arnol'd-Moser (KAM) barriers [1]. To demon- 
strate that this effect does not contribute to the observed localization in momentum space 
in Fig. 2 (a,c,d) of the main text (i.e. for atoms kicked by only a single optical lattice), we 
compare the experimental data and the simulated quantum trajectories of the per-particle 
energy ^ as a function of kick number A^ to simulated classical dynamics, averaged over 10^ 
classical trajectories. To briefly detail the determination of the averaged classical dynam- 
ics, each individual classical trajectory is determined by probabilistically projecting out the 
particle's momentum-space wavefunction after each application of the kick operator U. We 
show in FIG. A (i-iv) the classical dynamics of £:, as averaged over n = 10, 10^, 10^, and 
10^ individual trajectories, for the same system parameters as used in Fig. 2 (a) of the main 
text [T = 36 /xs ; r = 2 /xs ; Xi = 2.3 (^i = 50) ] K2 = 0]. For large sample sizes, the 
classical growth rate is found to be in good agreement with the expected classical diflFusion 
constant [2], As /AN - V^T^/2h^ = 2KI/k^, for the flrst lattice {2rfKl/K^ for the second 
lattice) . We have verifled that the slight deviation between the growth rate of simulated clas- 
sical trajectories and the expected classical diffusion is due to the flnite length (r = 2 /xs) of 
the lattice pulses. 




FIG. A. Averaged classical dynamics of the energy 5 as a function of kick number N for different 
sample sizes consisting of n trajectories, for off-resonant kicking with the first lattice [Ki = 2.3 
(51 = 50) ; i^2 = 0]. (i-iv) Averaged trajectories are shown for the cases of n = 10, 10^, 10^, 
and 10^ as solid red lines. In all plots, we also show as a dashed black line the expected classical 
diffusion for delta-function kicks. 



FLOQUET EIGENSTATE ANALYSIS OF U 

To gain a deeper insight into the expected response of the matter waves to a given kick 
operator U (i.e. for given pulse lengths, lattice depths, and pulse period), we analyze the 
properties of the matrix (7 in a time-independent manner. While monitoring the response of 
the matter waves as a function of the kick number A^ makes a direct connection to experi- 
ment and can discern between regimes of localization and delocalization in momentum-space 
through the long-term growth dynamics, at short times it may be subject to transient be- 
havior and small-scale fluctuations. We are experimentally restricted to a modest number of 
kicks to maintain a near-fleld treatment (assuming identical, overlapped spatial wavefunc- 
tions for all the matter- wave flelds) and to minimize contributions from nonlinear atom-atom 
interactions. Motivated by pioneering work that made a stroboscopic connection between 
the quantum 5-kicked rotor model and the ID Anderson model on a lattice [1], we study the 
properties of the Floquet (or Bloch-Floquet) quasi-energy eigenstates of the kick operator 
U. For regimes in which the growth dynamics result in dynamical localization, these Floquet 
eigenstates should all be localized in momentum space, while they are delocalized in the case 
of diflFusive growth. 

We flrst investigate the case of kicking oflF-resonantly with two lattices of equal strength, 
using the same parameter values as in Fig. 2 (e) of the main text, with variable stochasticity 
parameter Ki = K2 = i^i,2- In Fig- B (i,ii), we plot the momentum-space distribution of 
the lowest energy Floquet eigenstates of [/, 1^0^), for the values K12 = 2.15 and 2.34. The 
plotted distributions are shifted by the mean momentum p = (ipQlPlipQ) to p' = p — p. While 
the distribution is localized for the weaker kicking strength, it is delocalized into two regions 
for the larger strength. In fact, in Fig. B (iii) we observe a bifurcation of the eigenstate 
distributions as the stochasticity parameter K12 is increased. We can readily calculate the 
rms momentum width cr^/ of these distributions, and as shown in Fig. B (iv) there is a 
sharp rise in the width across a value of K12 ^2.2. Not surprisingly, this value is in general 
agreement with the observed onset of non-zero energy growth versus kick number in Fig. 2 (e) 
of the main text. 

Further motivated by the connection of our experiment involving atoms driven with two 
spatial frequencies to experiments performed with multiple temporal frequencies, which result 



in an effectively higher-dimensional localization when the frequencies are incommensurate, 
we can also perform a Floquet eigenstate analysis in the full parameter space of Ki and i^2- 
In Fig. C (i), we plot as white the regions in which the Floquet eigenstates are localized 
in momentum space and as dark blue the regions in which they are delocalized, with the 
criterion for delocalization being that more than 10% of the population resides further than 
5.5 momentum units (hki) away from the most populated mode. This plot shows localization 
along either axis, where dynamical localization is expected for any value of the kicking 
strength when only a single lattice is used, while delocalization generally occurs beyond 
some line in the parameter space spanned by Ki and i^2- If the expected connection to the 
scenario of driving with multiple temporal frequencies holds, this would be analogous to the 
phase diagram of the anisotropic 2D Anderson model, in which the pseudo-randomness of 
the site energies is provided by the fact that the pulse period does not match the Talbot time 
of either lattice, such that the relative kick strengths Ki and K2 give rise to off-diagonal 
terms that cause "tunneling" between momentum states. We also show in Fig. C (ii) the 
simulated change in e between one and A^ = 40 kicks, as in Fig. 3 (b) of the main text but 
over a larger range of Ki values. This plot demonstrates that qualitatively similar regions 
of localization and delocalization are also observed in the energy growth dynamics. 
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FIG. B. Time- independent analysis of the kick operator U for off-resonant kicking (T = 36 jus) by 
both lattices, with equal kicking strengths K12 for each. (i,ii) The momentum distribution of the 
lowest energy Floquet eigenstate of {7, \iPq)^ in the plane- wave basis. The momenta p are shifted 
by the mean value p = {tpQlPltpQ) to p' — p — p. For low values of the kicking strength, as in (i), 
the eigenstates are localized in momentum-space, while for larger values they became delocalized 
as in (ii). (iii) Momentum distributions as in (i,ii) as a function of the stochasticity parameter 
Ki^2' Beyond a value of i^i,2 ^ 2.2 a bifurcation and delocalization of the momentum distributions 
is observed, (iv) For the same range of stochasticity parameter values as in (iii), we plot the rms 
momentum width Gpi of the lowest energy eigenstates of {7, which exhibits a sharp increase across 
a value K12 ^ 2.2. 




FIG. C. Localization-delocalization transition in the K1-K2 plane, (i) As a function of the stochas- 
ticity parameters Ki and K2, we plot as white the regions in which the Floquet eigenstates of U 
are localized in momentum-space and as blue the regions in which they are delocalized, as detailed 
in the text. A general trend of delocalization for strong kicking with both lattices, and localization 
when kicking with only a single lattice, is observed, (ii) Here we replot the simulation data from 
Fig. 3 (b) of the main text, showing the change in energy e from the first to A^ = 40 kicks, but over 
a larger range of Ki values. Behavior qualitatively similar to the localization-delocalization plot in 
(i) can be observed in this dynamical response data. 
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